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Abstract 

We introduce new method of optimization for finding free parameters 
of afSne iterated function systems (IFS), which are used for fractal ap- 
proximation. We provide the comparison of effectiveness of fractal and 
quadratic types of approximation, which are based on a similar optimiza- 
tion scheme, on the various types of data: polynomial function, DNA 
primary sequence, price graph and graph of random walking. 



1 Introduction 

It is well known that approximation is a crucial method for making compli- 
cated data easier to describe and operate. In many cases we have to deal with 
irregular forms, which can't be approximate with desired precision. Fractal ap- 
proximation become a suitable tool for that purpose. Ideas for interpolation 
and approximation with the help of fractals appeared in works of M. Barnsley 
[2] and was developed by P. Massopust [6] and C. Bandt and A. Kravchenko [1]. 

Today we can apply fractals to approximate such interesting and interdis- 
ciplinary data as graphs of DNA primary sequences of different species and 
interbeat heart intervals [7], price waves and many others. 

Section [5] of this work is devoted to the construction of fractal interpolation 
functions. Necessary condition on free parameters di of affine iterated function 
systems is shown. One graphical example is given. 

In section[3]we give the common scheme of approximation of general function 
g e L'^[a, b] and obtain the equation for direct calculation of free parameters di. 

In section!?] we illustrate the results on concrete examples. 



2 Fractal Interpolation Functions 

There are two methods for constructing fractal interpolation functions. In 1986 
M. Barnsley 2; defined such functions, as attractors of some specific iterated 
function systems. In this work we use common approach, which was developed 
by P. Massopust [6]. 

Let [a, &] C R be a nonempty interval, 1 < TV G N and {{xi, yi) e [a, 6] x R | 
a = Xq < xi < ■ ■ ■ < xn-i < xn = b} — are points of interpolation. For all 
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i — 1,N consider affine transformations of the plane 



A■:M2^M^ A. 



(1) 



y ) \ J \ y J \ 

We require following two conditions hold true for all i: 

In this case 

_ Xj - Xj-i _ Ui - yi-i - di{yN - yo) 

— a — a 

_ bxj-i- axj _ by^^i - ayj - di(byo - ay^) 

~ u 1 Ji — I ; 

— a — a 

there {rfijili like family of parameters. Notice, that for all i operator Ai 
takes the line segment between {xo,yo) and {xN,yN) to the line segment passes 
through points of interpolation {xi_i,yi_i) and {xi,yi). 

Let /C be a space of nonempty compact subsets with HausdorfF metric. 
Define the Hutchinson operator |5j 

N 
i=l 

It is easily seen [2], that the Hutchinson operator $ take a graph of any con- 
tinuous function on a segment [a, 6] to a graph of a continuous function on the 
same segment. Thus, $ can be treated as operator on the space of continuous 
functions C[a, &]. 

For alH = 1 , TV denote 

Ui\[a,b]^[xi^i,Xi\, Ui{x) := a^x + Ci, , , 

Pi : [a, b] M, pi{x) := CiX + fi. 

Massopust [6] has shown, that $ acts on C[a, b] according to the rule 

N 

($g)(x) =^((p,ouri)(x)+d,(goM7i)(x))x[x,_i,x,](a;). (3) 
1=1 



Moreover, if |c?i| < 1 for all i = 1,7V, then operator $ is contractive on the 
Banach space {C[a, b], \\ ||oo) with contractive constant d < max{|di| | i = 1, N}. 
By the fixed-point theorem there exists unique function g* G C[a, 6], such that 
^g* — g* and for all g G C[a, b] we have 

hm 3*11^=0. 

n— >oo 

We will call g* fractal interpolation function. It is clear, that if g £ C[a,b], 
g{xo) = yo and g^x^) = y^r, then $((?) passes through points of interpolation. 
In this case we will call $"(5) pre-fractal interpolation functions of order n. 

Example 1 Picture shows fractal interpolation function, which was constructed 
on points of interpolation (0, 0), (0.5, 0.5) (1, 0) with parameters di = d2 = 0.5. 
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Figure 1: Fractal interpolation function. 

3 Approximation 



From now on we assume, that \di\ < 1 for all i = 1, TV. We try to approximate 
function g G C[a, b] by the fractal interpolation function g*, which is constructed 
on points of interpolation {(a^i, ?/i)}i^o- Thus, it is sufScient to fit parameters 
di G (—1, 1) to minimize the distance between g and g* . 

We use methods that have been developed for fractal image compression [5] . 
Notice, that from (0) and ([J) follows, that for all g,he L'^[a, b] 



W^g - ^h\\2 ^ \l ($3 - <i>/i)2 dx = 



N 



df / {g ou-^{x) - hou^^{x))^ dx 



Xi-i 



N 



< max||dj|| ■ . y^ai jg - h)'^ dx = max||d,|} ■ \\g - h\\2 

Thus, $ : L'^[a, b] — >■ i^[a, b] is contractive operator with a fixed point g*. 

Furthermore, instead of minimization of ||g — 5*||2 we will minimize \\g — 
^g\\2, that makes the problem of optimization much easier. The collage theorem 
provides validity of such approach. 

Theorem 1 Let (X, d) be a non-empty complete metric space. Let T : X ^ X 
be a contraction mapping on X with contractivity factor c < 1. Then for all 

xeX 

1 — c 

where x* is the fixed point of T . 
► For all integer n we have 

X*) < d{x, T{x)) + d{T{x), T^{x)) + • • ■ + d(r"-i(a;), r"(a;)) + d{T"{x), x*) 
< d{x, T{x)){l + c + c2 + . ■ . + c"-i) + d(r"(x), X*). 
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Letting n — > oo we establish the formula. 
Considering ([T]) and ([2]), we rewrite ([3]): 



i=l 

{yi - Vi-i)x + (xiyi-i - Xi-iyi) 



where 



ai{x) 



^ ^ iVN - yo)x + (a:,yo - x.^i^Af) 
(6 - a)a; + (xja - Xi„i6) 

= — 



Xi Xi—\ 

Thus, we have to minimize functional 

AT 



hf = X! / ~ "^(^) + di{Pt{x) - go 7j(x))) 



Ax. 



Setting partial derivatives with respect to di to zero we obtain 
("1(2;) - 9{x)) {Pt{x) - go 7,(x)) dx 



z-l,...,iV. (6) 



4 Discretization and results 

In this section we will approximate discrete data Z — {(z^, Wm)}m=ii a = ^^i < 
Z2 < • • • < zm = & by the fractal interpolation function g*, which is constructed 
on points of interpolation X = {(a;^, y^li^O' ^ ^ Taking X C Z, (xo,yo) = 
(zi, wi) and (a;jv,?/Ar) = (^;j\/,u'm) we fit parameters e ( — 1, 1) to minimize 

M 

^{Wjn - g*{Zm)f. 
m— 1 

Let us approximate Z by the piecewise constant function g : [a, b] — R. 
More precisely g{z) = Wm, where (z™, w„-i) G Z and Zm is a nearest neighbor of 
z. From ([6|) we obtain the discrete formulas for d^: 

- Wm) {f3i{z,n) - go 7i(z„i)) 

di — 2 7 z = 1, . . . , A* - 1. 

I] (^(^^m) -5°7i(^m)) 

Z7Ti£[2:i,a;i+i] 

(7) 

After finding d,; we obtain formulas for affine transformations Ai and we are 
able to construct fractal interpolation function g* for g. 
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Our aim is to compare fractal approximation with a piecewise quadratic 
approximation function which is based on the same discretization. On each 
segment [xi-i, Xi] approximating function has the quadratic form qi{x) — kix'^ + 
TiX + li. To get a continuous function we claim that qi{xi^i) = g{xi^i) and 
Qi{xi) — g{xi). From this we find coefficients ki and li. To find free parameter 
ri we minimize functional 

with respect to on each segment [xi-i, Xi],i = l,iV. The approximating 
function q{x) will have following form: 

{qi{x) — kix^ + rix + li, x e [a — xo,xi]; 

q2{x) = k2x'^ + r2X + I2, xe[xi,X2]; 
■ 

gAr(x) = kNx'^ + r^x + In, x € [xn-i,xn = b]. 

Since there is one free parameter in each function qi{x) and one parameter 
di for each afRne transformation Ai it makes the comparison correct. 

To compare fractal and quadratic approximations we consider four types of 
data. 

1. Polynomial function. 

2. DNA sequence. 

3. Price graph. 

4. Random walking graph. 

For all types of data M — 10000, — m, [a,b] — [1,A/], {'Wm}m=i are 
normalized sequences, that is E{{wm}) — and E{{w'^}) — 1. For all cases we 
choose {xo,yo) = {xN,yN) — {M,wm) and other interpolation points 

{xi, Ui), i = 1, N — 1 are local extremums of the given data. 

Example 2 Let f{x) = -6x + 5x'^ + 5x^ - 5x^ + x'^, x e [-1,2.5]. As we 
work with the segment [1,M] we map [—1,2.5] to it. Consider sequence Vm = 

f ^ 2(M-i) ^ ' "^^1' ^-^ ■ ''^m — {vm — Si)/s2, where si and S2 (iTe mean 
and deviation of {vm]^^^. Figure \^ shows the normalized sequence {wm}. 
Choose five interpolation points x^ — 1, xi — 500, X2 = 4000, X3 = 7500, 
Xi = 10000. Applying ^ we obtain di = 0.066, ^2 = 0.155, = 0.033, 
di = 0.096. The small values of \di\ mean that on segments [xi-i,Xi] fractal 
approximation function looks as a .straight line. Figure shows the graphs of 
fractal and quadratic approximating functions. 
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Figure 2: The graph of original function g. 




Figure 3: Fractal and quadratic interpolations of the polynomial function. 



Example 3 A DNA sequence can be identified with a word over an alphabet 
Af = {A,C,G,T}. Here we have the sequence of 10000 nucleotides of Edward- 
siella tarda. The graph represented by the formula 

P f +1, if m}'^'' nucleotide belongs to (A,G); 

' ^"^ 1 —1, if m^'^ nucleotide belongs to (C,T). 

For full description of representation of DNA primary sequences see '41- Figure 
\^shows the sequence {wm} after normalization of {vm} according to the formula 
in the previous example. Interpolation points are xq = 1, xi — 1000, X2 — 2500, 
X3 = 3000, = 3500,2:5 = 5000, X6 = 6500, xj = 7000, xs = 8000, = 9000, 
xiQ = 10000. Applying (Qj we obtain di = -0.001, ^2 = 0.274, = 0.31, 
d4 ^ 0.24, dg = -0.057, de = 0.211, dr = -0.42, ds = -0.121, dg = 0.215, 
dio — 0.158. Figure\^ shows the graphs of fractal and quadratic approximating 
functions. 

Example 4 We take price wave of 10000 prices v„n m — 1,M of one day 
period for EUR/USD, then normalize it (Figure\B^. Interpolation points are 
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Figure 4: Picture shows DNA Graph of 10000 nucleotides of Edwardsiella tarda. 




xo = 0, xi = 500, X2 = 1500, X3 = 2000, X4 = 2500,X5 = 3000,a:6 = 4000, 
X7 = 5000, = 6000, a;9 = 8000, xw = 10000. Applying we o6tom 
= -0.334, d2 = -0.004, da = 0.315, = 0.307, = 0.333, = -0.28, 
= -0.067, ds = 0.027, dg = 0.047, dio = -0.33. Figure\^ shows the graphs 
of fractal and quadratic approximating functions. 




Figure 7: Fractal and quadratic interpolations of the Price Graph. 



Example 5 Picture shows Random Walking Graph. It represented by the for- 
mula vq — 0, Vi = where is a random value with normal distribution. 
Interpolation points are xq = 0, xi = 1500, X2 = 2000, X3 = 3000, x^ = 4000, 



100 



50 




Figure 8: Normalized Random Walking graph. 

X5 = 5500, xe = 6300, X7 = 7600, xg = 8000, xg = 9000, xio = 10000. 
Applying (Qj we obtain di = -0.237, ^2 = 0.14, = -0.020, ^4 = -0.105, 
ds = 0.105, de = 0.0545, dr = -0.184, dg = -0.368, dg = 0.081, dio = -0.111. 
Figure shows the graphs of fractal and quadratic approximating functions. 

To compare the results we calculate approximation errors for each type of 
data. Let h{x) be the approximating function for data {wrn}m=i- Then ap- 
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Figure 9: Fractal and quadratic interpolations of the Random Walking graph. 



proximation error is 




{h{x,n) ~ W„i))2 



Here we represent the table of approximation errors for each type 



Fractal Quadratic 

PolynomialFunction 0.0359037 0.0245094 

DN A Primary Sequence 0.0692072 0.0624714 

Price Graph 0.0501345 0.0533686 

Random Walking 0.1015339 0.101438 



From it we see, that fractal approximation is better for price graph and nearly 
equal for random walking, but much worse for smooth function and slightly for 
DNA sequence. Different results were appearing during calculations of errors. 
We assume that some conditions could give us more exact approximation results 
from fractal interpolation function and for that extra observations should be 
established. 
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